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ABSTRACT 

We present Atacama Large Millimeter Array (ALMA) detections of atomic carbon line and dust 
continuum emission in two UV-luminous galaxies at redshift 6. The far-infrared (FIR) luminosities 
of these galaxies are substantially lower than similar starbursts at later cosmic epochs, indicating an 
evolution in the dust properties with redshift, in agreement with the evolution seen in ultraviolet (UV) 
attenuation by dust. The [C ii] to FIR ratios are found to be higher than at low redshift showing that 
[C il] should be readily detectable by ALMA within the reionization epoch. One of the two galaxies 
shows a complex merger nature with the less massive component dominating the UV emission and 
the more massive component dominating the FIR line and continuum. Using the interstellar atomic 
carbon line to derive the systemic redshifts we investigate the velocity of Lya emission emerging from 
high-z galaxies. In contrast to previous work, we find no evidence for decreasing Lya velocity shifts 
at high-redshift. We observe an increase in velocity shifts from z « 2 to z « 6, consistent with the 
effects of increased IGM absorption. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
high-redshift 


1. INTRODUCTION 

Observations from the Hubble Space Telescope have 
given us a broad-brush picture of the evolution of galax¬ 
ies over cosmic time. Galaxies in the reionization epoch 
at redshifts z > 6 were rarer, smaller, bluer and had 


the more evolved universe 1 

Ono et al. 

2013 

Bouwens 

et al.|2014b|a 

Smit et al.|20] 

1^. 'These observations can 


ogy with initial generations of stars forming out of low 
metallicity gas in relatively low mass halos. However, 
many details are still to be resolved. 

Whether the early galaxy population is able to gener¬ 
ate enough photons to complete cosmic reionization by 
redshift 6 is still a matter of debate and depends upon 
uncertainties su ch as the dumpiness of the intergalactic 
medium (IGM) ( Finlator et al.|2012|, the escape fra ction 


of ionizing photons from galaxies (Siana et al. 2015|_a2d 


the UV spec tra of very low-mass galaxies (|Schaerer & de 
Barros|2010 ). Strong evolution in the fraction o f galaxies 


showing Lyg emission (Pentericci et al. 2011 Schenker 


et al. 2012 


__ __ suggests a rapid change in 

the IGM neutral hydrogen fraction at z ~ 7. Some of 
the strength of this evolution could also be accounted 
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for by increased Lya abs orption within galaxies due to 
lower Lya velocity shifts (Choudhury et al. 20141. The 
first measurements of Lya velocity shifts in two z > 6 
galaxies do ind eed appear to show smaller shifts than at 
lower redshift (Stark et al.||2014a). 

Up til now most ol our knowledge of high-redshift 
galaxies is based on rest-frame ultraviolet (UV) and opti¬ 
cal observations. Observations in the far-infrared (FIR) 
can provide information on the physical conditions of the 
interstellar medium (ISM) and dust and molecular gas 
in obscured star-forming regions, as well as gas kinemat¬ 
ics. At z ^ 2 FIR photometry of rest-frame UV-selected 
galaxies shows that the star formation rate (SFR) is un- 
deresti mated by a factor of 5 from UV measurements 
alone (Reddy et al. 2012). With the incredible sensitiv¬ 
ity of the Atacama Large Millimeter Array (ALMA) it 
should be possible to detect dust continuum emission, 
molecular and atomic gas in high-redshift galaxies to 
complete our understan ding of their physical properties 
(Garilli & Walter 2013). In addition, the line data pro- 
vides information on the gas kinematics and can be used 
to derive dynamical masses, separate spatially unresolved 
mergers and define a systemic redshift against which to 
measure the emergent Lya kinematics. 

Several studies using ALMA and other facilities have 
resulted in non-detections and claims that the FIR line 
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on UV SFRs (Walter et al. 

2012 

Kanekar et al. 2013 

Ouchi et al.|2U13l|Ota et al. 

120141 |Gonzalez-L6pez et al.| 

2014 Maiolino et ai.||2015 

S 

chaerer et al.||2015). ALMA 


the fine-structure line of singly-ionized car bon, [Cii], in 


the z > 6.5 galaxies Himiko and IOK-1 (Ouchi et ah, 
20131 lota et al. 120141), despite the strong UV continuum 

■ III I I I I I < n i ^ ^ ^ ^ 1 1 


indicating Sh'R'^ 20 — lOOMQyr”^. They showed that 
the FIR (cool dust) contributions to the spectral energy 
distributions are similar to nearby dwarf irregulars, more 
than an order of magnitude below the levels expected 
for nearby starbursts/spirals such as M82 and M51. The 
physical interpretation of the lack of dust and [C ii] emis¬ 
sion in these z > 6 galaxies is a very low metallicity ISM. 
Observations of local low metallicity star-forming galax¬ 
ies such as I Zw 18 also show very low dust-to-gas and 
dust-to-stars ratios, further confirming such ga laxies as 
analogues of high-z galaxies (Fisher et al.||2014). 


This picture has brightened recently with seve ral detec¬ 
tions of the ISM in distant galaxies. Although [Maiolino 


et al. (2015) failed to detect three 6.8 < z < 7.1 Lyman 
Break Galaxies (LBGs) with ALMA in line or contin¬ 
uum, they did detect [C il] emission at the expected red- 
shift of BDF3299, offset by 4kpc from the UV position, 
which t hey interpret as an accreting or satellite clump 


of gas. Watson et al. (2015) detected a gravitationally- 
lensed galaxy at z = 7.5 in dust continuum from which 
they infer SFR-pm about 3 times greater than SFRjjy, 
similar to the ratio typically observed at z ~ 2 (Reddy 


et al.|20i2). The lack of a [C ii] detection in this galaxy is 


likely due to an uncertain redshift, as Watson et al. note 
the lack of an emission line redshift means the ALMA 
data only cover 50% of potential [Cii] line frequencies. 


Capak et al. (2015) observed nine LBGs at 5 < z < 6 
with ALMA and detected [C ii] emission from all of them 
and dust continuum in four galaxies. They showed that 
the dust emission in these galaxies is much weaker than 
expected based on similar galaxies at lower redshift. The 
Capak LBGs were selected based on interstellar UV lines 
so likely have higher metallicity than most previously- 
targeted high-z galaxies, although still lower than at low 
redshift. 

We are carrying out an ALMA program targeting UV- 
luminous LBGs at z > 6. In this paper we present obser¬ 
vations of the first two galaxies observed from our sam¬ 
ple. The galaxies were first identified in Canada-France- 
Hawaii Telescope (CFHT) optical and near-IR imaging 
and both have spectroscopic redshifts based on a con¬ 
tinuum break plus Lya emission. Th e galaxy which we 


name here as CLM 1 was discovered by Cuby et al. (2003) 
at a redshift of z = 6.17. The other galaxy, WM H 5 at 


z = 6.07, was discovered in Willott et al. (2013b). Both 
galaxies have near-lR continuum magnitudes AB ~ 24 
making them among the brightest rest-frame UV galax¬ 
ies known at this epoch. 

In Section 2 we describe our observations and in Sec¬ 
tion 3 our results. Section 4 discusses the ratio of [Cii] 
to FIR continuum and Section 5 the Lya velocity shifts 
observed. We draw conclusions in Section 6. Cosmologi¬ 
cal parameters of Ho = 67.8 km s~^ Mpc~^, Um = 0.308 


and Ua = 0.692 (Planck Collaboration et al. 2015) are 
assumed throughout. 


2. OBSERVATIONS 
2.1. ALMA 

Observations were made during ALMA cycle 2, in June 
2014, with between 29 and 32 antennas, and a maximum 
baseline of 650m. A total bandwidth of 7.2 GHz was 
employed in Band 6 using 4 dual-polarization sub-bands 
between 249 and 272 GHz. One of the sub-bands was 
centered at the [Cii] line (rest frame 1900.5369 GHz) for 
each source. The rest of the bands were used for con¬ 
tinuum measurements. A total of 95 minutes on-source 
integration time was obtained for each galaxy. 

The initial data editing and calibration were performed 
as part of standard data processing by the ALMA staff. 
The calibrated visibility data were then re-analyzed, per¬ 
forming additional flagging of bad time periods and bad 
channels. The data were re-imaged using the CASA 
Briggs weighting of the visibility data with Robust = 1 to 
create continuum and channel images. For both sources, 
the Gaussian restoring beam was close to circular, with 
a FWHM = 0.50”. Spectral cubes at 15.625 MHz resolu¬ 
tion (17.5 km s~^) were synthesized, and smoothed spec¬ 
trally for subsequent analysis, as required. The contin¬ 
uum was subtracted in the image-plane, using the off-line 
channels in the line cube. 

The spectra were transformed from the observed 
topocentric system to the local standard of rest (LSR). 
The velocity offsets applied were « — llkms”^. 

Spectral and spatial Gaussian fitting was performed 
using the CASA viewer, CASA fitting tools and custom 
software. Results for total line fluxes, continuum flux 
densities, and parametric source sizes, are given below. 
Flux uncertainties have 10% added in quadrature for ab¬ 
solute flux calibration uncertainty. 


2.2. Near- and Mid-Infrared 

To compare with the ALMA data we use near-i nfrared 
(NIR) data from the ESC VISTA VIDEO survey (Jarvis 
et al. 2013). Images from Data Release 3 in bands 
Z, Y, J, H and Ks were obtained, sampling the rest- 
frame ultraviolet continua. Astrometric calibration of 
the VIDEO images to the radio reference frame used 
by ALMA was performed by m atching bright VI DEO 
sources to the AllWISE catalog (|Cutri et al.||2013|) and 
extracting the ALMA phase calibrator positions from the 
AllWISE catalog. Given the size of the residuals this pro¬ 
cess gives an uncertainty on astrometric frame matching 
of « 071, comparable to the positional uncertainties of 
the NIR data due to S/N and seeing. Fluxes in the five 
NIR bands of the two target galaxies were determined 
using aperture photometry with aperture corrections to 
total fluxes. Due to the similar shapes, fluxes and S/N 
of the data in the three bluest VIDEO bands, a deeper 
combined ZY J image of each field was generated. Nei¬ 
ther galaxy is significantly spatially resolved in the NIR 
data. 

The two galaxies have been observed with the Spitzer 
Space Telescope in the I RAC 3.6 and 4.5/xm b ands dur¬ 
ing the SERVS survey (|Mauduit et al.||2012|). Due to 
the proximity of bright galaxies along the iine-of-sight 
to both z > 6 galaxies it was necessary to perform si¬ 
multaneous fitting of all objects in the IRAC field using 
the ZY J VIDEO image as a prior for the object s hapes. 
Models were fit to the NIR images using GALFIT (Peng 
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rest-trame U V continuum and [u ii] emission are co-spatial. i'wo toreground galaxies are visible on the right side ot the image. 


Observed Frequency (GHz) 
265.5 265.0 



Fig. 2.— ALMA [Cii] (upper) and VLT FORS2 Lya (lower) 
spectra of CLM 1. The ALMA spectrum is well fit by a single 
Gaussian emission line, with a marginal red wing excess. The Lya 
spectrum shows typical blue-absorbed asymmetry. Vertical blue 
and red dotted lines mark the [Cii] Gaussian fit peak and Lya 
flux-weighted centroid, respectively. Lya is offset from [C ii] by 
+430 ±69 kms-T 


et al. ]Mot and subs equent fitting of the IR AQ data was 
done with PyGFIT (Mancone et al. 2013). Both z > 6 
galaxies are clearly detected at both 3.6 and 4.5/rm, al¬ 
though with flux uncertainties of ^ 30%. 

2.3. Lyman-a spectroscopy 

Optical spectroscopy, including the detection of asym¬ 
metric Lya emission lines and continuum breaks, o f both 
galaxies were presented in their d iscovery papers (Cuby 
et al.||2003| |Willott et al.||2013b|). In order to compare 
these spectra to the ALMA spectra they were corrected 
from observed wavelength to the LSR by applying veloc¬ 
ity corrections of -1-26kms~^ and -1-32 kms~^ for CLM 1 
and WMH 5, respectively. The spectrum of CLM 1 was 
obtained using FORS2 at the ESO VLT with a spectral 
resolution of i? = 1400 and that of WMH 5 using CMOS 
at Gemini-North with R = 1000. 

3. RESULTS 
3.1. CLM 1 

For both galaxies there are strong detections of the 
[Cii] line evident in the datacubes at close to the ex¬ 
pected position and velocity. In Figure we present the 
ALMA and NIR imaging of the galaxy CLM 1. The 
[C ii] map is made from a sum over all channels showing 
signihcant line emission. The [Cil] emission is spatially 
extended along a similar direction as the beam exten¬ 
sion. A Gaussian fit to the source gives an observed size 
of Oi'TO X 0'.'42 at position angle east of north (PA) = 
42, compared with a beam size of 0"52 x 0"44 at PA = 
55. The CASA fitting software does not provide a de¬ 
convolved source size due to the only slightly larger size 
than the beam and outputs that the intrinsic source may 
be as large as (/.'Sd x C/.'OO. 

The white contours in Figure show the 1.2 mm con¬ 
tinuum emission as measured from the 3 sub-bands with¬ 
out the [C ii] line. In contrast to the strong and clear 
[C ii] emission the continuum in this source is very weak. 
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There is a marginal continuum detection at the 2 (t level 
that is only considered as plausible true emission be¬ 
cause its centroid co-incides to within (X'S of the [C ii] 
line centroid. There are several other peaks of this mag¬ 
nitude within the 6 " x 6 " image of Figure [l] that may 
be due solely to noise. We note that there is an emerg¬ 
ing trend in ALMA observations of non-ULIRG high- 
redshift galaxies that the [Cil] line emission is detected 
at much higher significance than the continuum, despite 


the large instantaneous bandwidth of 7.5 GHz (Riech- 
ers et al.]|2014[ [Willott et al.||201^ 2015[ Gapak et al. 


20151. We discuss the ratio of line to continuum emission 
further in Section [d] 

The NIR image in Figure shows the location of the 
rest-frame UV continuum in GLM 1. The Lya contri¬ 
bution to this flux is negligible. The galaxy is indistin¬ 
guishable from a point source at the O'.'S resolution of the 
NIR data. The centroids of the [C ii] and rest-frame UV 
co-incide to within Cf.'l, suggesting they trace the same 
star forming regions of the galaxy, at this resolution. 

Visible to the east of GLM 1 is a lower redshift ellipti¬ 
cal galaxy whose centroid is only 3" from GLM 1. The 
presence of this galaxy combined with the high UV lumi¬ 
nosity of GLM 1 raises the possibility that the high lu¬ 
minosity is in part due to magnification by gravitational 
lensing. To determine the potential lensing magnification 
we have analyzed the properties of this elliptical galaxy. 
We fit galaxy models to 12 band photometry from the 
GFHT Lega cy Survey, VIDEO and AllWISE using the 


(jkriek et al. 

12009 

ED) is well 

ht by 


"By an old stellar population 
at z = 0.54 with a stellar mass of 5 x 10^° Using 

the observed correlation betwe en velocity dispersion and 
stellar mass (Wake et al.||2012| at z = 0. 1 and with neg- 


ligible evolution from z = 0.5 to z = 0.1 ( Shu et al.|2012 1 


this gives a velocity dispersion of 126 kms“^. Assuming 
an isothermal sphere potential this lensing configuration 
provides a magnification factor of only 1.13. Given this 
relatively low magnification we do not make any lensing 
corrections to physical values in this paper. 

Figure plots the [C ii] spectrum of GLM 1 and the 
Lya spectrum for comparison. Interpretation of the ve¬ 
locity difference between [G ii] and Lya is deferred to Sec¬ 
tion!^ The [G ii] spectrum is well fit by a single Gaussian 
plus a very low level, flat continuum. The fit continuum 
level is consistent with the low S/N continuum detec¬ 
tion from the other 3 sub-bands described above. There 
is marginal excess emission in the red wing. From the 
Gaussian peak frequency we determine the systemic red- 
shift of the galaxy as z = 6.1657 ± 0.0003. The Gaussian 
FWHM is 162 ± 23 kms“^. Table 1 contains further 
measurements made from the ALMA data. 

With this velocity and the upper limit on the source 
size (< 0"55 or 3.2 kpc) we can calculate an uppe r limit 
to the dyna mical mass. Following the procedure of jWangj 



Fig. 3.— Observed-frame optical to far-IR SED of CLM 1 (blue 
circles). The best-fit model from GIG ALE is shown as a black 
curve and the model fluxes through the observed filters as red cir¬ 
cles. The grey curves are 39 nearby galaxy SEDs with UV to far-IR 
photometry from the KINGFISH survey ([Dale et al.||2012|| rang¬ 
ing from dusty infra-red galaxies to dwarf irregulars with very low 
thermal dust emission. The lower panel shows the residuals from 
the best-fit model. CLM 1 has a SED most similar to the least 
dusty nearby dwarf irregulars. 


TABLE 1 

Millimeter data for the z > 6 LBGs 



CLM 1 

WMH 5 


6.176 ±0.002 

6.076 ±0.001 

2[CII] 

6.1657 ± 0.0003 

6.0695 ± 0.0003=" 
6.0639 ±0.0003*" 

FWHMfcii] 

162 ± 23 kms-i 

251 ± 15 kms-l =" 

68 ±13 kms-l*" 

l[Cii] (Jykms-i) 

0.234 ±0.031 

0.659 ±0.071 

-^[Cii] (h©) 

(2.40 ± 0.32) X 10® 

(6.60 ±0.72) X 10® 

/l.2mm (p-Jy) 

44 ±26 

218 ±41 

hpiR (h©) 

(2.5 ± 1.5) X 10i° 

(1.25 ±0.24) X 10** 

h[CII]/hpiR 

(9.5 ± 5.7) X 10-® 

(5.3 ± 1.2) X 10“® 

SFR[cn, (Moyr-l) 

24 ±3 

66 ±7 

SFRFiR(AL©yr-i) 

7±4 

33 ±6 

SFRsed (M©yr-l) 

37 ±4 

43 ±5 


Notes.— 

The redshift and FWHM of component ‘A’ only. 

The redshift and FWHM of component ‘B’ only. 

Uncertainties in Lpip and SFR only include measurement uncer¬ 
tainties, not the uncertainties in extrapolating from a monochro¬ 
matic to integrated luminosity, luminosity-SFR calibrations, or un¬ 
certainty in star formation history. 

optical attenuation and its re-emission in the infrared. 
Stel lar lig ht is modeled with a single jBruzual &: GharT] 
lot (2003) stellar population with constant star forma¬ 


et al. 


tion rate, Ghabrier initial mass function and metallicity 
= 0.008. Dust attenuation followed the ICalzetti et al.l 


1 IS tne unknown inclination angle. 

IR via the SED parameterization of 

Dale et al. 

(2014 

)■ 

Using photometry from the rest-frame U V, optical and 

'I'he highest likelihood parameters are a ste 

lar mass 


far-IR we ht galaxy spectral models to understand the 
evolutionary state of GLM 1. A synthetic ALMA contin¬ 
uum filter is generated using the wavelengths covered by 
the three line-free sub-bands. We us e the Python imple- 
mentation of the CIGALE package ( Roehlly et al.|[20T2 1 
to fit the SED. This code determines the ultraviolet- 


of 1.3 X 10^^ Mq, SFR= 37M 0 yr“^ and minimal dust 
attenuation of E(B-V)=0.0I. This very low amount of 
dust attenuation is constrained by the blue rest-frame 
UV slope of /3 = —2.0 (where/A oc A^) and in an IR 
remission scenario consistent with the low FIR luminos¬ 
ity from the marginal 1.2 mm continuum detection. We 
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note that an increasing star formation rate and/or higher 
fraction of the IRAC flux contributed by nebular emis¬ 
sion lines could substantially decrease the stellar mass, 
bringing it more into line with the dynamical mass upper 
limit. It is also possible that the [C ii]-emitting gas is not 
tracing the full dynamical mass of the system. 

In Figure ^we show the observed photometry of CLM 
I along with the best-fit model photometry, the model 
spectrum and, in the lower panel, the residuals from 
the fit. With almost as many free parameters as data 
points, a good fit has been found. The observed IRAC 
fluxes are slightly higher than the model suggesting possi¬ 
bly stronger nebular lines than modelled, although there 
are large uncertainties in the IRAC fluxes. The gray 
curves in Figure 1^ show redshifted SEDs for nearby 
galaxies from the HNCFISH survey, normalized at rest- 
frame 150 nm. These galaxies have been observed with 
GALEX, Spitzer Space Telescope, Hers chel Space Obser- 


SFR in very low dust and/or metallicity galaxie s (Ouchij 


et al. 2013 Fisher et al. 2014 Ota et al. 2014). An al- 


vator y and in the optical and near-IR (Dale et al. 2007 
20121 . The galaxies range in total IR luminosity froni 
10^ to 10^^ Lq. CLM 1 has a SED most similar to the 
^ 10^ Lq galaxies, despite having an IR luminosity of 
> 10^° Lq, in essence it is a scaled up (by 10^) version 
of the low metallicity dwarf irregulars. 

It has recently been shown that the [C ii] luminosity 
is an effective tracer of the star formation rate in low 


redshift starbursts (De Looze et al 


2014|). In addition to the far-IR luminosity and SED 


l.|2014[ 
L lumin 


Sargsyan et al. 


(largely constrained by UV, rather than FIR, fluxes) this 
gives three independent measures of the SFR in this z > 
6 galaxy. For the [C ii] luminosity we use the relation 
SFR (MQyr“^) = 1.0 x 10“^L[cn] (Aq) (Sargsyan et al. 


2014). 

We calculate the far-IR luminosity Lfir (integrated 
over rest-frame 42.5 to 122.5 pm) from the observed 
1.2 mm continuum assuming a greybody spectrum with 
dust temperature, Td = 30K and emissivity index, (3 = 
1.6. We use Td = 30K because that is the dust temper- 
ature for similar FIR-lu minosity galaxies at low-redshift 
(Symeonidis et al. 20131. We note that a higher dust tem- 
perature of 7 d = 45 K would increase the Lfir value by 
a factor of about three. To convert from Lrir to SFR we 
use the relation SFR (Mq yr“^) = 1.5 x 10~^°Lfir, (Lq ) 
appropriate for a Chabrier IMF (Carilli & Walter||2013 ). 

At these redshifts the Cosmic Microwave Background 
(CMB) has a temperature of 19 K and can potentially 
bias measurements of dust continuum luminosity. There 
are two competing effects: (i) a high CMB background 
against which the continuum is measured, and (ii) an 
increase in the dust temperature due to heating by the 
CMB. Since we only have one continuum point and no 
constraints on dust temperature, we cannot make an ac¬ 
curate correction for t hese effects, but n o te tha t accord¬ 
ing to the analysis of |da Cunha et al. (2013) the two 
effects are of comparable size and opposite sign for likely 
dust temperatures at this redshift, so we make no cor¬ 
rection. 

All three SFR estimates are listed in Table 1. For CLM 
1 we find that SFRser and SFRjqh] are comparable, but 
SFRfir is lower than SFR[cn] (see Section and lower 
than SFRsed, as expected from the low dust contribu¬ 
tion to the SED in Figure|3 Previous studies have shown 
that SFRfir may be unreliable as a tracer of the total 


ternative is that a higher dust temperature of 1 d = 45 K 
would raise SFRfir to a level comparable with SFRsed 
and SFR[cii]. Future observations at shorter wavelength 
are critical to constrain the full IR SED of high-z galaxies 
and constrain the dust temperature. 


3.2. WMH 5 

The ALMA and NIR images for WMH 5 are shown 
in Figure]^ There is a much more firm 1.2 mm con¬ 
tinuum detection in this source at significance 6cr. The 
beam is more circular for WMH 5 than for CLM 1 and 
both the [C ii] and 1.2 mm continuum are clearly spa¬ 
tially extended at PAw 80. The [Cii] emission centroid 
shows a significant offset from the NIR emission of Of'4. 
The NIR emission is consistent with a point source. The 
Lyo emission is also spatially unresolved along the slit 
direction. 

Figure shows the [C il] spectrum of WMH 5. The 
line can be split into two Gaussian components. We la¬ 
bel the component with higher flux and larger linewidth 
(FWHM=251 kms“^) as ‘A’ and the other component 
as ‘B’ (FWHM=68kms“^). Assuming the linewidths 
trace mass gravitationally we identify the systemic red¬ 
shift of the main galaxy WMH 5 with the velocity of 
‘A’. To determine the nature of this multiple veloc¬ 
ity component system we show in Figure a position- 
velocity diagram along the major axis of [C ii] emis¬ 
sion. The separation of the two components in velocity 
is clear, but we also find a spatial offset of (y.'3 (2kp(A 
between the peaks of the two components. In Figure 
we mark the centroids of components ‘A’ and ‘B’ sep¬ 
arately to highlight how those positions compare with 
each other and the NIR emission. Both [Cil] compo¬ 
nents are marginally resolved with intrinsic major axis 
sizes of 0"48 ± 0"1 and (/.M ± 0"1 for ‘A’ and ‘B’, respec¬ 
tively. We determine dynamical masses from the [C il] 
sizes and FWHM of M^yn = 1-7 x 10^°/sin^iMQ for A 
and Mdyn = 1-0 x 10®/sin^fMQ for B. 

We carry out SED-fitting using CIGALE for the entire 
WMH 5 system. The observed and modelled SEDs are 
shown in Figure]^ As for CLM 1, the photometry is well 
fit by a single stellar population whose attenuated UV 
emission is re-radiated in the FIR. The highest likelihood 
parameters are a stellar mass of 2.3 x 10^° Mq, SFR = 
43MQyr“^ and dust attenuation of E(B-V)=0.05. The 
higher dust attenuation in WMH 5 is constrained by the 
redder UV spectral slope and consistent with the higher 
Lfir. However, in this system we know that there is a 
spatial offset between the UV and FIR, so such a simple 
scenario of attenuation and re-emission is not physically 
plausible. 

From this analysis the system appears to be an on¬ 
going merger of two galaxies. ‘B’ is apparently lower 
mass but is spatially coincident with the NIR emission. 
This is similar to som e lower redshift ULIRG mergers 
(Chapman et al. 2004) with an optically obscured com- 
ponent that dominates the FIR emission and a lower 
Lfir component less obscured by dust. Similar cases 
of most of the [C ii] associated with optically-faint com- 
ponen ts have been reported at z > 5 by [Capak et al.| 
(2015); Maiolino et al. (2015), however in our case com- 
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Fig. 4.— Left: The background image is ALMA integrated [Cii] line map of WMH 5. White contours show the clear 1.2 mm dust 
continuum detection (contour levels 2, 3, 4, 5 cr bearn”^). The NIR centroid (black plus symbol) is significantly offset to the east. 

Right: The background is the 2 :y J NIR image. Contours show the ALMA [Cii] emission from the left panel at contour levels 3, 6, 9 cr bearn”^. 
Grey crosses labeled ‘A’ and ‘B’ correspond to the centroids of the two [C ii] velocity components. The rest-frame UV continuum is consistent 
with ‘B’, the component with lower [Cii] luminosity and velocity width. There is a foreground galaxy further east of the system. 


Observed Frequency (GHz) 

269.5 269.0 268.5 



Fig. 5.— ALMA [Co] (uppeU and Gemini CMOS Lya (lower) 
spectra of WMH 5 (see Figure!^ for details). The [Co] data have 
been continuum-subtracted in ^e image plane. For this galaxy 
there are two [Cll] velocity components (‘A’ and ‘B’), that are 
each well fit by a Gaussian. The systemic redshift is taken to be 
the Gaussian peak of ‘A’ as this is the most luminous and highest 
width component. The Lyo spectrum is asymmetric and the flux- 
weighted centroid offset by -1-265 it 52 kms“^ from ‘A’ and -1-504it 
52 kms“^ from ‘B’. 



Fig. 6.— Position-velocity map for the [Cll] emission in WMH 
5. The NIR galaxy centroid is located approximately along the 
major axis of the [C ll] emission and its position is identified by the 
solid black line labelled NIR. The two [Gll] components, ‘A’ and 
‘B’ are clearly separated in position and velocity, with ‘B’ close to 
the NIR centroid, as seen in Figure]^ 


ponent ‘A’ also dominates the FIR continuum and has a 
large linewidth so is not some peripheral, satellite gas 
cloud. The stellar mass estimate of ^ 2 x 10^° Mq 
is > 5 times larger than the dynamical mass of ‘B’ of 
1.0 X 10®/ sin® i Mq, unless the inclination angle i < 30°. 
Either ‘B’ is a disk viewed extremely face-on or it is not 
responsible for all the stellar mass. 

For WMH 5 the SFR estimates from the far-IR, [C ii] 
and SED-fitting in Table 1 are all within a factor of two. 
This is rather surprising given the complex nature of the 
system. The millimeter and UV SFR have similar values, 
but their spatial displacement indicates that they are not 
probing the same star-forming regions, so the true SFR 
should be the sum of both components. 
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Fig. 7. — Observed-frame optical to far-IR SED of WMH 5 (blue 
circles), best-fit CIGALE model and SEDs of nearby KINGFISH 
survey galaxies (see Figure[^for details). WMH 5 has a SED most 
similar to nearby dwarf irregulars. 

4. THE [Cii] - FAR-IR LUMINOSITY RELATION 

For both [C ii] and FIR luminosities to act as reli¬ 
able star formation rate indicators requires that the ra¬ 
tio of the two behave in a predictable manner without 
too high a scatter. At low-redshift most galaxies have 
ratios in the range 10“^ to 10“^ with a fairly linear re- 


2014 

Sargsyan et al. 

to 

o 

the ULiRG regime w 

here a 


2014). The exception to this is in 


usually observed, thought to be related to extreme densi- 
ties and temperatures in circumnuclear starburst regions 


dFarrah et al.|20f3 Mngdis et al.|2014 |Gonzalez- Alfonso, 

et al.||2015p. Figure]^ plots the low redshift data with 

small black circles. 

At higher redshift (z > 4), the first [Cil] measure¬ 
ments were made in very high Lfir sources and showed 
low [Gii]/FIR ratios comparable wit h nearby ULIRGs 
(Maiolino et al.|2005 lono et al.|2006 ). Subsequent work 
at redshifts between 1 and 7 reyealed a wide range of 
ratios from 10“"* to 10“^. At the low end the re are high- 
luminosity z > 5.7 quasars (Wang et al. 2013) and hyper- 


luminous infrared galaxies (HyLlRGS; Lfir > 10^^ Lq) 


(Riechers et al.||20i3 2014). At the opposite end of the 
range, high [GiiJ/hlR ratios up to 10“^ have been found 
in some massive 1 > z > 2 ULIRG s lying on the main- 


sequence of star-fo rming galaxies (Stacey et al. 2010_ 
Brisbin et al.|201^ ). In these galaxies, the star-formation 


is more spatially extended than in nearby ULIRGs lead¬ 
ing to higher [C Ii] luminosities. 

In Figure we also include recent ALMA observa- 
tions of Lftr objects comprising z > 6 AGN 

dWillott et al.||2013a[ |2015|) and z > 5 galaxies ([Capak 
et al. |2015[ |Maiolino et al.|2015| plus GLM 1 and WMH 
R Our data are the hrst non-limit data at z > 6 in 
-^FiR < 10^^ Lq galaxies without AGN. For z > 5 galax¬ 
ies the [Gii]/FIR ratios display a broad range at some¬ 
what higher values than at low-redshift. Ratios of ^ 10“^ 
suggest extended star-formation with low metallicity and 
an intense radiation field. 

The similarity of high-redshift galaxies hosting quasars 



Fig. 8. — Ratio of [C 
versus far-IR luminosit}' 
are plotted as red square 
in the [C ii] transition a 
to-right) BDF3299 Clur 
and HFLS 3 at 2 : = 6 
5 < 2 < 6 fmagenta circ 
et aL| ||2014||; |Capak et £ 
luminosities have been d 
hosting galaxies at 2 : > 

i] to far-IR (42.5 to 122.5/im) luminosity 
. The two 2 > 6 LBGs from this paper 
s with error bars. Other galaxies detected 
b 2 > 6 blotted as orange circles are (left- 
Qb A at 2 = 7.107 (iMaiolino et al.||2015i 
337 (IRiechers et al.||!i(Jl!J|). (galaxies at 

es) are trom|Rawle et al.| ||2Ui4|); |Riechers 
Ll.||j2015|l where the (Japak et al. total IK 
viaed bv a factor 1.75 to give Lt?tr . AGN- 
5.7 are shown with blue stars dMaiolino 

et aL|2005l IVenemans et al. 120121 IWang et al. 120131 |\Villott et al. 

il013al Galaxies 

(green circles) are Irom 

(including some AGJN) at 0.5 < 2 < 5 
Brisbin et al.|(|2015|l and the compilation 

ofIDe Looze et al. (12014 

' and at 2 < 0.4 (black plus svmbol) from 

the compilation ol |Gracia-(Jarpio et al.| (2011 and in prep.), ihe 
2 > 6 LBGs CLM 1 and WIVIH b have ratios comparable to 2 > 5 


galaxies of similar Lpir, somewhat higher than at 2 < 0.4. 


and without quasars gives us confidence that the star 
formation properties are only weakly affected by black 
hole accretion onto moderate mass black holes of Mbh 
1 O®A70, compared to the very low [Gii]/FIR values 
found in the more luminous (in bot h accretion and FIR ) 
z > 5 Mbh ~ 10® Mq quasars of Wang et al. ([2013). 
We note that the Lfir < 10^® Lq z > 6 AGN have ob- 
served ratios /i.2mm/.l[cii] similar to GLM 1 and WMH 
5, their lower ratios in this figure being due to the factor 
of 4 higher Lfir calculated using a dust temperature of 
Td = 47 K, compared to = 30K for the galaxies. 


5. Lya EMISSION LINE VELOCITY SHIFTS AND 
INTERPRETATION OF RAPID EVOLUTION IN Z ^ 7 
Lya GALAXIES 

Atomic or molecular gas in star forming regions spread 
throughout galaxies is ideal for measuring the systemic 
redshift of the galaxy. At high redshift such measure¬ 
ments have been used to determine the ionized bubble 


sizes surrounding luminous quasars (Garilli et al. 2010). 


In this Section we use this information to determine the 
velocity shifts of the observed Lya emission from GLM 1 
and WMH 5. Lya at z > 6 is the most important probe 
of the fraction of neutral hydrogen in the IGM and the 
process of cosmic reionization. 

In Figuresand 1^ we plotted the [G ii] and Lya spectra 
on the same velocity scale after correcting the datasets 
to the local standard of rest. For both galaxies the Lya 
line is asymmetric with a broader red wing than blue and 
the line center is shifted to the red from the systemic 
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redshift defined by the [Cil] observations. This highly 
asymme tric line shape is chara cteristic of Lya at high- 
redshift ( Shimasaku et al.|2006 ) due to neutral hydrogen 
absorption of the blue wing. For CLM 1 the measured 
offset is +430±69 kms“^ and the Lya flux drops to zero 
just before the [C ii] line center. For WMH 5 there is un¬ 
certainty in the systemic redshift due to the two compo¬ 
nents. If ‘A’, the component with larger [Cii] linewidth 
and therefore likely higher mass, is to be identified as the 
Lya systemic redshift then the offset is ±265±52 kms“^. 
However, the astrometry suggests that ‘B’ is the origin 
of the rest-frame UV and therefore Lya emission and in 
this case the shift is ±504±52 kms“^. For the remainder 
of this section we consider both cases as possible. 

The emergence of Lya from galaxies is a complicated 
process that involves resonant scattering off neutral hy¬ 
drogen and absorption by dust. The observed Lya pro¬ 
files depend upon factors such as geometry, gas covering 
factor, dust and outflow velocity. Irrespective of the ion¬ 
ization state of the IGM, galaxies tend to show Lya pro¬ 
files offset to the red due to resonant scattering and selec- 
ti ve absorption w ithin galactic-scale outflows (see review 


Dijkstra|2014 ). At 2 ; > 6 where the neutral fraction of 


the IGiVI becomes significant, it is expected that IGM ab¬ 
sorption of Lya decreases the over all Lya flux and shift s 
the line center further to the red (Laursen et al. 2011). 
The very sharp decrease in Lya line strength of LBGs 
between z = Q and z = 7 has been interpreted as evi¬ 
dence for a rapid change in IGM neutral fraction at this 
epoch, in tension with reionization models that predict a 
smoother change in neutral fraction over this short cos- 


Treu et al. 112013). 


mic time rtPentericci et al. 2011 Schenker et al. 2012 


The hrst measurement of Lya velocity shifts at red- 
shifts relevant to the reionization epoch have recently 
been made by|Stark et al.| (|2014a|). These authors used 
the UV nebular GiiiJA1909 doublet, which is a rela- 
tive ly strong line in you ng, low-mass star-forming galax¬ 
ies (Stark et al. 2014b), to determine the systemic red- 
shifts of two strong Lya emitters. The gravitationally- 
lensed z = 6.027 galaxy A383-5.2 h as a high Lya equiv- 
alent width with strong asymmetry. Stark et al. (2014a I 
quote the Lya peak velocity offset as 120 kms“^ 


t; 


or 

comparison with our velocity offsets and previous mea¬ 
surements at lower redshift we determine from their Fig¬ 
ure 5 an approximate centroid-based velocity offset for 
this galaxy of 150±30 kms“^. The z = 7.213 galaxy GN- 
108036 has a more moderate Lya equivalent width and a 
mu ch lower significance detection of Cm] was obtained 
by Stark et al. (2014a I with a very low and uncertain 
v elocity ottset of —37 ± 113 kms“^ 


Stark et al. (2014a) noted that these offsets are small 
compared to those of LBGs at lower redshift and this 
may be due to different physical conditions at high red¬ 
shift. In particular there is known to be a negative cor- 
relation between velocity shift and Ly a equivalent width 
( Shibuya et al.|20l4 Erl] et al.||2014 ) and Lya strength s 


increase with redshift up to z = 6 (Stark et al. 2011). 

Since Lya more easily escapes from galaxies via the wings 
than the line center, a decrease in velocity shift at high- 
redshift could be partially responsible for the sudden de¬ 
crease in Lya emitters at this epoch. [Choudhury et ah 
(2014) showed that a negative evolution of velocity shift 
with redshift can match the observed Lya equivalent 



A383-5.2 


2 < 2 < 3 .i Galaxies (Erb et al.] 
t^z>6 Galaxies (Stark et al.) 
^ 2>6 Galaxies (this paper) 
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Fig. 9.— Rest-frame equivalent width of the Lya emission line 
versus the velocity offset between Lya and the systemic redshift. 
The two z > 6 LBGs from this paper are plotted as red squares 
with error bars. For WMH 5 the two [Cii] components, ‘A’ and 
‘B’, provide two estimates of the velocity offset connected by a 
dotted line. Magenta circles show the only other two galaxies at 
such h igh redshifts with measured Lya velocity offsets (|Stark et al.| 
|2014a|| based on the offset to the Clll]A1909 UV emission line. A 
compa rison sample of galaxies at intermediate redshift from |Erb| 
|et al.| | [20T4| also Strom et al. in prep.) is plotted as green circles, 
it the systemic redshift of the WMH 5 Lya emission is given by 
‘B’, then the two z > 6 galaxies from this paper lie at the large 
velocity shift end of the distribution defined by the intermediate 
redshift sample. 

width distribution evolution to smooth reionization his- 


tory models (see also Bolton & Haehnelt 2013). 

In Figure!^ we compare the Lya velocity shifts mea¬ 
sured for CLM 1 and WMH 5 to those of A383-5.2 and 
GN-108036. We choose to plot the shifts as a function 
of observed Lya rest-frame equivalent width because of 
the correlation between velocity shift and Lya equiv¬ 
alent width. In addition to these z > 6 galaxies we 
show for comparison a sample sel ected at lower r edshift 
(2 < z < 3.1) as LAEs and LBGs ( |Erb et al.|2014 Strom 
et al. in prep.). For the lower redshift sample, spectro- 
scopic equivalent widths are plotted when available, oth¬ 
erwise photometric equivalent widths are used. Several 
things are evident from looking at Figure]^ The shifts 
for CLM 1 and WMH 5 are larger than for A383-5.2 
and GN-108036. None of the high-redshift galaxies have 
shifts at odds with those observed at 2 < z < 3.1 when 
compared to galaxies with similar Lya equivalent width. 
If we adopt ‘B’ as the systemic redshift for WMH 5, 
then both galaxies in this paper have Lya velocity shifts 
at the high end of the distribution, consistent with ad¬ 
ditional IGM abs orption at highe r redsh ift as predicted 
by the models of Laursen et al. (2011). Our analysis 
shows no observational evidence for velocity shifts that 
decrease as a function of redshift, at a given Lya equiva¬ 
lent width. Obviously a much larger sample of observed 
high-redshift galaxies compared to realistic simulations 
incorporating radiative transfer would be required to use 
the Lya velocity shift distribution to constrain the de¬ 
tails of reionization. 

6. CONCLUSIONS 
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We have presented ALMA [Cii] line and dust contin¬ 
uum detections of two UV-luminous LBGs at redshift 
z > 6. These detections were made in relatively short 
integrations in Early Science operations providing hope 
that other galaxies in the reionization epoch with lower 
SFR will be detectable with the full ALMA array. Our 
results, in accord w ith other recent ALMA observations 
(Capak et al. 2015), confirm that, despite unexpectedly 
low Lpift, the increase in the [Cii]/FIR ratio at high- 
redshift leads to [G ii] lines that are bri ght and very use¬ 
ful tra cers of the ISM, as predicted by [Walter fc Garilli| 
( 2008|). The resolution of the merger of WMH 5 with 
a fairly compact ALMA array illustrates the power of 
[C il] line observations to understand the details of star 
formation and galaxy assembly at this epoch. 

CLM 1 and WMH 5 were selected for study because 
their UV luminosity, stellar mass and redshift are compa¬ 
rable to the galax y Himiko which was unexpectedly unde¬ 
tected by ALMA (Ouchi et al.|2013| Ota et al. 20141. The 
main difference between our target galaxies and Himiko 
(and most other high-z non-detections) are our relatively 
low Lya equivalent widths. Since the Lya line strength 
decreases with increasing dust (and hence metallicity), 
this could explain the difference between these results. 
We note that the recent success with ALMA at 5 < z < 6 
is also for UV-lumin ous galaxies with relatively weak Lya 
( Capak et al.|[20T5 ). 

However, the UV slopes of our galaxies are fairly blue, 
indicating little dust absorption of the UV photons and 
we do observe a deficit of FIR photons compared to ex¬ 
pectations, as shown by the SEDs that are similar to 
nearby dwarf irregulars, not with nearby galaxies with 
similar SFRs of tens of solar masses per year. Future 
observations at shorter wavelength are critical to con¬ 
strain the full IR SED of high-z galaxies and constrain 
the dust temperature. Are they very cool like low metal¬ 


licity dwarfs, despite the higher UV photon density and 
CMB temperature? Determining this is important as the 
full IR SED is required to derive more accurate values of 
SFRpin to determine whether this is a more or less ac¬ 
curate SFR indicator at high-z than SFRjcn] or SFRuv- 
We have for the first time used millimeter lines to de¬ 
termine Lya line velocity offsets from systemic for galax¬ 
ies in the reionization era. We find that the shifts in our 
two galaxies are comparable to the high end of the dis¬ 
tribution at lower redshift, consistent with predictions 
from Ly g line asymmetry du e to neutral hydrogen ab¬ 
sorption (Laursen et al. 2011). More detailed studies in 
the future of spatially-resolved kinematics of [C ii] and 
Lya in such galaxies will allow stronger constraints to be 
placed on the escape of Lya emission from galaxies and 
its subsequent absorption by the IGM. 


Thanks to staff at the North America ALMA Re¬ 
gional Center for processing the ALMA data, Jean- 
Gabriel Cuby for providing the optical spectrum of 
CLM 1, Dawn Erb for providing unpublished data on 
intermediate redshift Lyman Break Galaxies, and Pe¬ 
ter Capak and Caitlin Casey for interesting discussions. 
This paper makes use of the following ALMA data: 
ADS/JAO.ALMA#2013.1.00815.S. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF 
(USA) and NINS (Japan), together with NRG (Canada) 
and NSC and ASIAA (Taiwan), in cooperation with the 
Republic of Chile. The Joint ALMA Observatory is oper¬ 
ated by ESO, AUI/NRAO and NAOJ. The National Ra¬ 
dio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. 

Facility: ALMA. 
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